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Abstract. For positive p-harmonic functions on Riemannian mani- 
folds, we derive a gradient estimate and Harnack inequality with con- 
stants depending only on the lower bound of the Ricci curvature, the 
dimension n, p and the radius of the ball on which the function is de- 
fined. Our approach is based on a careful application of the Moser iter- 
ation technique and is different from Cheng- Yau's method 2 employed 
by Kostchwar and Ni [5], in which a gradient estimate for positive p- 
harmonic functions is derived under the assumption that the sectional 
curvature is bounded from below. 



1. Introduction 

The study of harmonic functions on Riemannian manifolds has been one of 
the central subjects in geometric analysis. In their classical work Cheng- Yau 
[15j derived the following gradient estimate for positive harmonic functions 
on Riemannian manifolds: 

Theorem A (Cheng-Yau) Let M be an n- dimensional complete Riemann- 
ian manifold with Ric > — (n — 1) k, where k > is a constant. Suppose 
that u is a positive harmonic function on a geodesic ball B (o, R) . Then 

(1-1) sup < C n , 

B(o,R/2) u K 
where C n is a constant depending only on n. 

An important feature of Cheng- Yau's estimate is that the RHS (which 
stands for the right hand side) of (II. ID depends only on n, k and R, it 
does not depend on the lower bound of the injectivity radius or a global 
coordinate system. From PDE viewpoints, deriving a Harnack inequality 
requires some bounds on the coefficients in some fixed coordinate system, 
thus not suitable for many problems defined on manifolds. We also observe 
that the RHS of (jl.ip is optimal in the sense that even for k > 0, the bound 
stays bounded when R — > oo. 

There are two major ingredients in the proof of Theorem A. First a 
Bochner formula is used to derive a lower bound of the Laplacian of | Vit| 2 for 
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a harmonic function u in terms of the lower bound of the Ricci tensor. The 
second major ingredient is a clever application of the maximum principle. 
The trick is to multiply |Vu| 2 by a cut-off function, derive a new differen- 
tial inequality for the product and then apply the maximum principle. The 
cut-off function is constructed using the distance function. As a result, the 
new differential inequality involves the Laplacian of the distance function. 
As is well known, the Riemannian distance function is uniformly Lipschitz 
and its Laplacian has an upper bound depending on the lower bound of the 
Ricci tensor. 

Cheng- Yau's approach turned out to be very useful and some important 
results for other problems are deeply influenced by Theorem A. For example, 
P. Li [8j obtained the sharp lower bound for the first eigenvalue of a manifold, 
which was later generalized by Li-Yau [9 J . Similar results were also obtained 
by Li-Yau [10] for heat equations. S. Y. Cheng [1] and H. I. Choi [3] obtained 
gradient estimates for harmonic mappings, etc. We refer to |15t HI] and the 
more recent survey [7] for an overview of the subject. 

p-harmonic functions are natural extensions of harmonic functions from 
a variational point of view. It has been extensively studied because of its 
various interesting features and applications. Compared with the theory for 
harmonic functions the study of p-harmonic functions is generally harder 
because the equation, even though elliptic, is degenerate and the regular- 
ity results are far weaker (see, for example [H]). Recently, there has been 
renewed interest in p-harmonic functions. In particular R. Moser [11] es- 
tablished a nice connection between p-harmonic functions and the inverse 
mean curvature flow. In a recent paper [5] Kotschwar and Ni derived, among 
other things, a local gradient estimate for p-harmonic functions under the 
assumption that the sectional curvature is bounded from below. It is re- 
markable that the constant in their estimate does not blowup when p — > 1, 
which leads to interesting results on the inverse mean curvature flow prob- 
lems. Their proof follows the same strategy introduced by Cheng- Yau [2] for 
harmoinc functions (i.e. p = 2). However for general p-harmonic functions, 
the computation involves the full Hessian of the distance function when the 
cut-off function is introduced. As a result, a lower bound on the sectional 
curvature has to be assumed in [5]. 

Kotschwar and Ni speculated that their estimate may hold if only a lower 
bound on the Ricci tensor is assumed. The main result of this paper is to 
establish the following theorem: 

Theorem 1.1. Let (M n ,g) be a complete Riemannian manifold with Ric > 
— (n — 1) k. Assume that v is a positive p-harmonic function on the ball 



B(o,R) C M. Then there exists a constant C PyTl such that 
< C p>n (1 + y/RR) /R on B(o, R/2). 
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The proof of Theorem 11.11 will be presented in section two. As far as 
the second major ingredient of Cheng- Yau's proof is concerned, our ap- 
proach follows a different strategy by carefully using the Moser iteration 
technique. This approach only involves differentiating the distance function 
once and hence bypasses the difficulty of handling the full Hessian of the 
distance function. In the special case p = 2, when p-harmonic function are 
just harmonic functions, Theorem II .11 is exactly Cheng- Yau's theorem. An 
immediate consequence of Theorem 11.11 is the following Harnack inequality. 

Theorem 1.2. Let (M n ,g) be a complete Riemannian manifold with Ric > 
— (n — 1) K. Assume that v is a positive p-harmonic function on the ball 

B{o,R) C M. Then there exists a constant C Ptn such that for any x,y £ 
B(o,R/2), 

v(x)/v (y) < e CvA l +^ R ). 

It follows that if Ric > 0, then we have a uniform constant c Pjn (indepen- 
dent of R) s.t. 

(1.2) sup v <c pn inf v. 

B(o,R/2) ' B(o,R/2) 

This was already proved by Rigoli, Salvatori, and Vignati [12j. In fact, 
they proved the stronger result that fjl .2 j) holds provided that the volume is 
doubling and a weak Poincare inequality holds. See also [3]. 

Another standard application of Theorem 11.11 is the following Liouville 
theorem, which was also deduced from the Harnack inequality in |12j . 

Let u be a p-harmonic function bounded from above or below on a com- 
plete Riemannian manifold with non-negative Ricci tensor, then u is con- 
stant. 

Finally we point out that our constant C p ^ n in Theorem 11.11 becomes 
unbounded as p — >■ 1, while in Kotschwar-Ni's result, all the constants stay 
bounded when p — > 1. We do not know if the method can be tweaked to 
remove this defect. 

2. The gradient estimate for the p-harmonic functions 

p-harmonic functions arise naturally as critical points of the LP (p > 1) 
norm of the gradient. Let (M n ,g) be a complete Riemannian manifold and 
SlcMan open set. A function v G W, *? (f2) is p-harmonic if 

div (JVu| p ~ 2 W) = 

in the weak sense, i.e. 

/ |Vv| p ~ 2 (Vt>,vo = o 

Ju 
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for all £ G W ' p (0). By [16] for example, v must be C 1 '". Moreover v G 
if j> > 2; u G if 1 < p < 2. Away from {Vi> = 0}, v is in fact smooth. 
Suppose that is positive. Set u = — (p — 1) logv. Then n satisfies 

(2.1) div(|Vu| p - 2 W) = |Vn| p . 

Let / = |Vn| 2 . We define 

£ ty) = div (r' 2 - l A (VV)) - pf p/2 - 1 (Vu, VV) , 

where 

, . , , n . Vn ® Vn 
^ = z d + (p - 2) 5—. 

|Vn| 

We need the following lemma from [5] and the proof is by direct calcula- 
tion. 



Lemma 2.1. 



£ (/) = 2/p/ 2 - 1 (\D\\ 2 + Ric (Vn, Vn)) + (f - l) | V/| 2 /^ 2 

Remark 2.1. Lemma \2. 1\ holds point-wisely in {x : f(x) > 0}. From the 
gradient estimate of [16] we know that f = |Vn| 2 G C a /or some a > ane? 
/ G /or some j3 > 1. 

We choose a local orthonormal frame {et} with ei = Vn/|Vn|. Then 
(|2.ip takes the following form 



(p - 1) nn + ^ njj = /. 

i=2 

Therefore 



D 2 u\ >u 2 n + 2J2< + J2 u ii 

i=2 i=2 

>n 2 1 + 2^n 2 4 + -l- K>J 

i=2 \i=2 / 

n . 

2 1 +2^n 2 i + ; ^- T (/-(p-l)nii) 2 

i=2 

1 2 2(p-l) , ^ , (P-1) 2 ^ 2 2 

— rl 7- Mi + 1 + r~ u n + 2 >^ u i 

n — 1 n — 1 \ n — 1 I ^-^ 



n 2 



«=2 



1 2 2(p-l) ^ 2 

> r/ -Mi+ao) % 

n — 1 n — 1 *— i 



2 

i' 

8=1 
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where ao = 1 + min ( , 1^ > 1. Using the identities 

2 fu u = (Vu,V/>, 



v„ 2 _ 1 i v /i 

Z^ u ii - i—f 

we end up with 
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\D\$ > -L-/» - ^ <V„, V/) + 

1 1 n — 1 re — 1 4/ 

Assume that Ric > — (re — 1) n. Therefore 

(2.2) £ (/) > -2 (re - 1) ^/ 2 + (*±^° - l) |V/| 2 Z^ 2 

+ ^-r/ p/2+1 - ^fcil/^-i ( V n, V/) 
n — 1 n — 1 

> -2 (n - 1) fcf / 2 + -^-f / 2+1 - fcil /^-i ( Vu , V/) 
re — 1 n — 1 

Equation (|2.2p holds wherever / is strictly positive. Let -ff = {x G : 
/(x) = 0}. Then for any nonnegative function ip with compact support in 
ft \ K, we have 

(2.3) j (p ,2 - l Vf +(P- 2)P' 2 ~ 2 (Vu, V/> Vu, V^) 

+P [ / p/2 - 1 (Vu,V/)^ + ^- / r/ 2+ v 
Jn n-1 J n 

< 2(n-l)* / / /^(V^V/)^. 

In particular, let e > and V = / e 6? ? 2 where / e = (/ — e) + , r/ G Cq°(Bji) 
is non-negative, 6 > 1 is to be determined later. Then direct computation 
yields 

VV> = bf^Vfr] 2 + 2/^Vr?. 



Using the above in fj2.3[) we obtain 

(2.4) bj^ (/ p/2 - 1 /e^ 1 |V/| 2 + (p-2)/^ 2 (Vu,V/) 2 / e '- 1 )r / 2 

+ 2(p - 2) / /f/ 2 - 2 (Vu, V/) /^(Vu, Vn) + 2 / f^fy (Vf, V V ) 
Jb r Jn 

+ pf / p/2 - 1 (Vn,V/)/ e V + ^- / f/ 2+ W 

n-lj n 

< 2(n-l)« / /^-W + / /^(V^V/)/* 

Jn n-1 J n 
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Since u 6 C 1 - , / e C a and V/ G L^tt) for some a > and > 1, 
we see that except for the first term, all the other terms converge to the 
corresponding form without e. For the first term, observe that 

/P /2-l / 6-l| V/ |2 + {p _ 2)fP /2-2 (VU) v/) 2 fh -l > ai/p/2 _ 1/ 6_l |v/|2 

where ai = 1 if p > 2 and a\ = (p — 1) if p € (1, 2). Thus by passing e to 
we have 

(2.5) 6a! / /^-2| V/ |V 

+ 2(p-2) f f E ^ +b (Vu,Vf} V (Vu,V7 1 } + 2 [ f^viVf^v) 
Jn Jn 

+ P [ f^ +b (Vu,Vf) V * + ^- [ f^+V 
Jn n - 1 J n 

< 2(n - l)k [ /^ + V + / f^ +b <V«, V/) r? 2 . 

Jn n-1 J n 

From now on we use ai, a 2 , • • • etc. to denote constants depending only 
on p and n. Combining terms in (|2.5j) using the definition of / we have 



(2.6) a x b [ fP/^\Vf\W + / f / 2+1 +V 

Jn n - 1 Jn 

< 2(n-l)«/" / p/2+ V + « 2 / /^^IV/I?? 2 
Jn Jn 

+a 3 / f^+^lV/HVr^. 
Jn 

For i?3 (the third term on the RHS, Li, L 2 , i?i etc. are understood similarly) 
in (|2,6p we have 

/^/ 2 - 2 |v/| V + x / |Vr?l 2 / 6+p/2 . 

& Jn 



IP l< aib 

\R 3 \ < — 

Also by Cauchy's inequality i? 2 of (|2.6p can be estimated as 

\R 2 \<^ [ /P/2+b-2| V/ | 2?? 2 + ^5 /• / 6 +P / 2+ l 7? 2_ 

4 Jn o Jn 

With these two inequalities we obtain 

(2.7) ^ / r/ 2+6 - 2 iv/iy + ^- / f/^+y 

2 Jn n-lj a 

< f ( 2 (n - w + ^ivy) r/ 2 + fe + £ / /^/ 2 +y . 

By requiring 
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we see that the last term in on the RHS of (12. 7h is majorized by the last 
term on the LHS. Therefore we have 



^ l / p/2 + 6-2 |v/| 2 2 + 1 f /P /2+l 

< 2(n - 1)« / /p/ 2 +V + ^ / |Vr?| 2 / p / 2+b . 

o Jn 

For the first term on the LHS we use 

| V(/P /4+6/ 2r?) |2 < l ( P +6) 2 / p/^2| V/ |2 T/ 2 + 2/ p/2+6| V , 7 |2 < 

From the above we obtain 
(2.9) 

f iVC/^+^^p+d! / fP/ 2 +l+ b v 2 < K d 2 [ / p / 2+ V+«7 / lVr/1 2 / 6 ^/ 2 . 
Jn Jn Jn Jn 

where d\ ~ 6, c?2 ~ b (recall b > 1, d\ ~ 6 means (ii is comparable to 6, 
c?2 ~ ^ is understood the same way). 

The following Sobolev embedding theorem of Saloff-Coste plays an im- 
portant role in our approach: 

Theorem B (Theorem 3.1 of |13j ) Let (M n ,g) be a complete Riemannian 
manifold with Ric > — (re — 1) k. For n > 2, there exists C , depending 
only on n, such that for all B C M of radius R and volume V we have for 
feC^(B) 

j l/l 29 ) 17 ^ e C'(^R) V -2/n R 2 (^j |V/| 2 + iT 2 / a ) S 

where q = n/(n — 2). For n = 2, the above inequality holds with n replaced 
by any fixed n' > 2. 

From now on, we assume f2 = B (o,R). Theorem B gives 
(2.10) / f^=^rj— 



where co(n,p) > depends only on n,p. Let &o = ci(n,p)(l + \/7lR) with 
c\{n,p) large enough to make bo satisfy (j2.8j) . then (j2.9j) and ()2.10j) combined 
gives 

(C „, „ \ (n—2)/n 

/ f^ 1 ^) +a 8 bR 2 e c ^V- 2 / n / /P/2+i+y 

Jn / Jn 

< a 9 b 2 be C2bo V- 2 / n [ fP/ 2+b v 2 + a w e C2ho V- 2 / n R 2 [ \Vrj\ 2 f p / 2+b 
Jn Jn 
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Lemma 2.2. Let b\ = (&o + |)n^2- Then there exists 03(71, p) > such that 
( 2 - 12 ) H/ll^(^ /4 )^ c 3|^. 

Proof of Lemma 12. 2t Let b = bo in (|2.11|) . then by comparing L2 and Ri 
of (]2.1ip we observe that 

a 9 ^/ p/2+6 ° < \a,boR 2 f' 2+l+b » 

if / > anb^R" 2 . Thus in the evaluation of R\ we decompose Q into two 
subregions, one over the places where / < anb^R^ 2 and the second region 
is the complement of the first region. With this decomposition we have 

Now (|2.1ip with b = bo can be written as 

1 /oil \ \ (n— 2)/n 



(2.13) Qf Z^^ 1 ^^) n + | 6o i? 2 e- b ^- 2 /" jf f/^+V 

C2 fe yl-2/n + aioeC2 b y-2/n jR 2 /" | Vt? | 2 /^/ 2+&0 

./n 



p+26o 



Now we choose rj to make i?2 in f|2.13j> dominated by the LHS. Let rji G 
Cq 00 ^) satisfy 

< 771 < 1, r/! = l in 5 3i?/4 , |Vr7i| <C(n)/R. 

Let r/ = 77™ where m = 60 + | + 1. Direct computation shows 

2b +p 

(2.14) i? 2 |Vr/| 2 < ai 3 6or/ b o+p/2+i. 

By (|2.14p and Young's inequality, the R2 of (|2.13p can be written as 

1 2 fbo+p/2 



a w R 2 j lVr/1 2 /" 

C 2b +p 

< a u bl / / 6o+p/ Vo+p/2+i 

< a u bl(jj b ^+w) 



bp+p/2 
b +p/2+l 



y60+P/ 2 +l 



, ,fe +p/2+2 
< ^0 R 2 fb0+P /2+l 2 b 0+P /2 bo 

- 2 ./j/ '15 #26 +p 

With the estimates on Ri , i?2 we arrive at 
(2.15) ( / f(bo+p/2)n/[n-2)\ < a^e^V^Hl ( % 
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Recall bi = Oo + f)^- Taking the l/(b +p) root on both sides of ([215} 
we have 



\LH {B3R/A) <a 17 V*bZ/& 
Lemma 12.21 is established. □ 



Now we go back to (|2.1ip . by ignoring L 2 we have 
(2.16) U f n - 2 rjw=2\ 

To apply the Moser iteration we set 

ti R R 

bl + 1 = bl ^2> 0/ = B(o, 2" + ^), « = 1, 2.. 

and choose r# G (f2) s.t. 

7]i = 1 in n l+x , rji = 0'mQ\ni, |V^| < — , < 7# < 1. 
Then in ()2. 16j) , by letting b + | = , 77 = 77; we have 

tL /t,+ ')^ £ (""w)*(/ n , ( * +tf|vi » |2)/ f • 

By the estimate of |V77;| 

/ e c 2 b \ 1 / 6 ' , \ 1/6; 

(2-17) < ("ir^J (6g6, + 16') \\n L H { n x) . 

Notice that 57 = then (gJ7j) leads to 

(2.18) ||/||l~( Sr/2 ) 

/ e c A \Er =1 iA»/ / n V V /bl 

* { a ^v^j n & o (^2) + 16« j 11/11^^ 



" C2 i '0 



- ai9 T7T7bT 6 o 2i>1 ll/lli 6 i(B 8 «/4)- 

Using Lemma 12.21 in (I2.18P we obtain 

(2-19) \\f\\L°°(B R/2 ) < a 2 ob 2 /R 2 . 

Thus Theorem 11.11 is established. □ 
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